Topological Structure of a Vortex in Fulde-Ferrell-Larkin-Ovchinnikov State 
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We find theoretically that the vortex core in the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state 
is quite different from the ordinary core by a simple topological reason. The intersection point 
of a vortex and nodal plane of the FFLO state empties the excess spins. This leads to observable 
consequences in the spatial structure of the spontaneous magnetization. We analyze this topological 
structure based on the low lying excitation spectrum by solving microscopic Bogoliubov-de Gennes 
equation to clarify its physical origin. 

PACS numbers: 74.25.Op, 74.25.Jb, 74.70.Tx 



Much attention has been focused on exotic supercon- 
ductors of various kinds, such as heavy Fermion materi- 
als, high T c cuprates, or organic conductors. These sys- 
tems are characterized by strongly correlations and/or 
low-dimensional electronic structures. Here pair poten- 
tial changes the sign in the momentum space for the rel- 
ative coordinates of a Cooper pair, such as p-, d-, /-wave 
symmetries while s- wave pairing does not. 

The Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state 
Q differs from these varieties, since the pair poten- 
tial changes sign in the real space for the center of mass 
of a Cooper pair. Thus the superfluid phase in the FFLO 
state is described by a spatially modulated order parame- 
ter. The FFLO state is widely discussed as for its realiza- 
tion in various research fields, ranging from superconduc- 
tors in condensed matter , neutral Fermion superfluids 
in an atomic cloud Q , to color superconductivity in high 
energy physics [jj. Originally the FFLO state is proposed 
when a superconductor is under an applied field where 
up and down spin populations become unequal by the 
Zeeman energy shift 111, 12, l| , followed by several theoret- 
ical works 0>M El Ell- The excess unpaired carriers, say 
up spin electrons, are accommodated in the nodal posi- 
tion of the order parameter in space by forming a midgap 
bound state near the Fermi level |2|. 

A key factor to stabilize the FFLO state is the relative 
weight between the Pauli paramagnetic effect and orbital 
depairing effect due to screening current under a field, 
apart from the material aspect such as purity and so on. 
Usually the latter overwhelms the former, thus th e upper 
critical field H C 2 is orbitally limited. CeCoIns [ll|, LLJ, 
El 0, E3 is a prime candidate for the strong Pauli 
effect because of the first order phase transition at H C 2 
in higher fields. In fact, there are several pieces of strong 
evidenc e, p oin ting to the FFLO state in this material 

In order to better characterize the FFLO state and to 
help its identification, we undertake a detailed study of 
a vortex state whose topological structure turns out to 
be quite distinct from conventional flux line. This study 
is indispensable because the FFLO state may only be 
induced in vortex state in a superconductor. It is also 



pointed out that the associated quasi-particle structure 
around a vortex core is directly probed experimentally 
via scanning tunneling spectroscopy (STM-STS) through 
local density of states, and neutron diffraction or X-ray 
diffraction through the induced magnetization modula- 
tion. 

We consider here a single vortex in a superconductor 
along the z-axis with cylindrical symmetry and the order 
parameter changes its sign at the z = plane. We solve 
the Bogoliubov-de Gennes (BdG) equation with the Pauli 
paramagnetic effect through the Zeeman effect. The Zee- 
man effect plays a critical role to stabilize the FFLO mod- 
ulation plQ- We fully take into account both orbital and 
spin effects in an equal footing 0, ^| . This enables us 
to reveal an interesting topological structure in a FFLO 
vortex. 

We start with the BdG equation for the quasi-particle 
wave functions u q (r) and f q (r) labeled by the quantum 
number q: 
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with JC = — h 2 V an d the self-consistent equation 

A(r)=g ]T Wq (r)^(r)/( £q ), (2) 

where /(e q ) = \/{e £ ^ kBT -f 1) is the Fermi-distribution 
function, and g (< 0) is the coupling constant. We 
set the energy cutoff uod = 5Ao. To prepare the un- 
equal spin population for a =|,|, the chemical poten- 
tial is shifted as fi^i = \i ± Sfi where the plus (mi- 
nus) sign corresponds to the up- (down-) spin. We 
assume the axially symmetric system with a radius R 
(= 2£o) and a height 2Z (= 5£o)> and isotropic s- 
wave symmetric pairing. Hence, we write the eigenfunc- 
tions as i£ q (r) = ?x q (r, z) exp [i(qe — \)0] and v q (r) = 
v CL (r : z)exp[i(q e + \)0] with A(r) = A(r, z)e~ w for a 
single vortex in cylindrical coordinates, qe is an an- 
gular momentum of the eigenstates. Along the radial 
and axial directions, the rigid boundary conditions are 
imposed: u q (r = R,0,z) = v q (r = R,0,z) = and 
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FIG. 1: Spatial structures of Re[A(r) / Ao] (a) and m(r, z) / po 
(b). po is the density in the uniform BCS state at T = 0. Note 
that m(r,z) is missing near the origin r = z = 0. 

w q (r,6>,z = ±Z) = v^(r,0,z = ±Z) = 0. Throughout 
this paper, A /e F = 0.12, <5/i/A = 0.28, and T = are 
fixed, where Ao is the energy gap in the BCS state at 
T = 0. We performed full self-consistent computations 
under these conditions. 

The spatial dependence of the order parameter 
Re[A(r)] is shown in Fig. 1(a), where the vortex line 
threads at r = and the nodal plane of the FFLO state 
is at z = 0, as shown schematically in Fig. 2. A(r) van- 
ishes on the vortex line at r = and on the nodal plane 
at z = where the sign of A(r, z) changes. It is seen that 
A(r) approximately varies sinusoidally as a function of z 
in this example. 

In Fig. 1(b) we show the spatial structure of the spin 
polarization m(r) = (r) — p;(r), which is defined with 
the density of each spin component p a as 

m(r) = £ [S(r)| 2 /(e q ) - K(r)| 2 (l + /(e q ))] • (3) 
q 

It is seen that m(r, z) is induced both along the vortex 
core at r — and on the nodal plane at z = because 
the order parameter vanishes there. However, the spin 
polarization is absent near the origin (r = 0, z = 0) where 
the FFLO nodal plane (z = 0) and the vortex line (r = 0) 
intersect. This can be understood in terms of the phase 
structure of the order parameter A(r) around r = and 
z = 0. 

We display schematically the relative configuration of 
the vortex line and nodal plane in Fig. 2. Quasi-particles 



tracing either the path (i) or (ii) feel the sign change of 
A(r) by cutting through the vortex line with 27r-phase 
winding (i) or by crossing the nodal plane (ii) (see also 
Fig. 1(a)), giving rise a bound state near the Fermi level 
because the order parameter's phase is shifted by tt on 
these paths. This is the so-called 7r-shift physics 0. 
The former (latter) is the core-localized (node-localized) 
state. These states accommodate excess up-spin elec- 
trons and yield spontaneously the spin polarization at 
these particular locations. In contrast, the quasi-particles 
tracing the path (in), which goes through the origin, do 
not experience the sign change because of tt + tt phase 
shift, one coming from the vortex contribution and the 
other from the nodal contribution. Thus it yields no low- 
energy state due to 7r-shift, meaning the absence of the 
spin polarization at the origin. Thus the spatial struc- 
ture of the spin polarization shown in Fig. 1(b) is deeply 
rooted to the topology of the FFLO vortex. Note in 
fact that these characteristics are absent for the vortex 
with the winding number 2 (Air- vortex) by obvious reason 
(27T + 7T shift). 

We now consider the electronic structure associated 
with the FFLO vortex, namely the local density of states 
(LDOS), N(j(r,E) which is calculated by 

N^(r,E) = J2K(r)\ 2 5(E-e cl ), 
q 

N l (r,E) = J2Mr)\ 2 S(E + e (i ). (4) 
q 

In Fig. 3, we show the spectral evolution N^(r : E) 
along the paths A(r = 0, z = 0) -> B(0,Z/2) -> 
C(R/2,Z/2) D(R/2,0) A presented in Fig. 2. We 
plot the spectrum N^(r,E) and N^(r^E) at the posi- 
tion A, B, and D in Fig. 4. The energy gap opens in 
\E — < Ao in N^(r,E). Outside the gap, the con- 
tinuum spectra are seen. At the vortex core outside the 
FFLO node (Fig. 4(a)), N } (r,E) and A/j(r, E) have, re- 
spectively, the zero-energy peak at /i| and /ij due to the 
core-localized state. In the path B — > C in Fig. 3, we 
see a typical spectral evolution found in the conventional 
vortex . As we move away from the core site B to- 
wards C, the higher angular momentum state with the 
higher energy acquires the spectral weight. On the other 





\ \ 


^ B \ 




M 






— r 

i ' A 


Id i c t 




— > 



FIG. 2: Schematic view of the vortex line and nodal plane. 
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hand, in the nodal plane of the FFLO state, the peak 
of the core states for up (down) spins at ^ is split 
into two. Along the vortex line A — > B in Fig. 3, we see 
that the split core-localized bound states tend to merge, 
forming a usual core bound state at B. Away from the 
core site A towards D along the nodal plane, the splitting 
of the two peaks decreases. In the FFLO node without 
vortices, these two peaks are merged at E = /xj . 

The LDOS structures N^(r,E) and N^r.E) are 
closely related to the local spin polarization m(r) in 
Fig. 1(b). At the vortex core (Fig. 4(a)) or at the nodal 
plane (Fig. 4(c)), most of the peak states for up (down) 
spins are located at E < (E > 0), and they are occupied 
(empty). Thus, the spin polarization appears. However, 
at the vortex core in the nodal plane (Fig. 4(b)), since 
the one of the split peaks is occupied both for up and 
down spins, the up and down spin populations are equal. 
This is the reason from the electronic structure why the 
spin population is absent in the vortex core on the nodal 
plane. 

Obviously STM-STS is a direct method to see these 




FIG. 3: Spectral evolutions of iVf (r, E) along the paths (A — > 
B — > C — > D — > A) shown in Fig. 2. Density plot (upper) 
and stereographic view (lower) are displayed. In the energy 
scale, E — corresponds to the Fermi energy /i. 
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FIG. 4: The local density of states for spin up (solid lines) 
and spin down (dashed lines) components at the positions B 
(a), A (b), and D (c). 



spectral evolutions in the FFLO state. In particular, 
comparison of tunneling conductance at A and D shown 
in Fig. 4 yields direct evidence for the FFLO state. 

The origin of the splitting of the vortex core bound 
state at the position A is due to the Doppler shift by 
the FFLO modulation. If we assume A(r) = Aoe ±l( ^ zZ 
to simply catch the essential feature, the wave functions 
can be written as u q (r) = w q (r, z )e ±t( ^ zZ ^ 2 and w q (r) = 
v q (r, z)e ±zQzZ / 2 which satisfy Eqs. (1) and (2). Thus, in 
the kinetic term K in Eq. (1), V 2 Z -> V 2 Z ±iQ z V z - Q 2 z /4,. 
Since the quasi-particles traveling across the nodal plane 
with the finite z-component vp, z °f the Fermi velocity are 
important to form the bound states at the nodal plane, 
the Doppler shift term ±iQ z V ' z ex Q z vf, z gives the up- 
per or lower energy shift of the spectrum, producing the 
split of the peaks for the vortex core states. We also cal- 
culate the LDOS spectrum using the test potential with 
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A(r, z) = A(r) sin (Q z z). It is seen that with increasing 
the splitting of the peak becomes wider. 

We should comment on the NMR experiment of 
CeCoIns where "novel" high field phase for H || c is iden- 
tified as the FFLO state. Our result is basically consis- 
tent with the data: Double peaks in the resonance spec- 
trum are related to the spin polarization structure. One 
peak is associated with the nodal plane of the FFLO 
state, the other is for the bulk far from the nodal plane. 
Detailed calculations show that we obtain a double-peak 
shape in the magnetization distribution profile, corre- 
sponding to the Knight shift in the NMR experiment. 

We propose here neutron diffraction experiment on the 
FFLO state. If we look at a scattering plane (0,q y ,q z ) 
of the reciprocal space (H,K,L) for External || z, there 
is no magnetic Bragg diffractions at (0,0, ±Q Z ) because 
the scattering vectors are parallel to the spin polariza- 
tion (|| H). However, around the higher Bragg points 
at (0, 27rn/L, ±Q Z ) (L is the vortex lattice spacing, n = 
1, 2, 3, • • •) we have the spin polarization component per- 
pendicular to the scattering vectors, which can be observ- 
able in this particular setting for diffraction experiments. 
Note that the vortex lattice spacing L ~ 200A for IT and 
Q~ x is a long period (Q^ 1 changes widely depending on 
H and T). That must be small angle scatterings. The 
scattering images the spin polarization, consisting of the 
sheets with the period of Q z along the field direction and 
rod-like objects along vortex lines with voids of the po- 
larization at the intersection. This spatial modulation 
of the spin polarization can be detected by synchrotron 
X-ray diffraction via magneto-striction, where the lattice 
modulation is caused by the coupling to the spin modu- 
lation in the FFLO nodal plane. 

As mentioned before, the present results are indepen- 
dent of the detailed material parameters because the 
topological structure in the FFLO vortex is essential. As 
for the pairing symmetry; We assume here isotropic s- 
wave pairing. The d-wave pairing which is advocated 
for CeCoIns 18] does not alter our conclusions based 
on the 7r-shift of the FFLO state. In that case, the 
extra-nodal excitations associated with vortex core at the 
Fermi energy extend towards the nodal directions of the 
anisotropic pairing function [l^ . 

We remark that the present considerations can be ap- 
plicable to resonance Fermionic superfluid under rotation 
with mismatched Fermi surfaces of two species in ultra- 
cold atoms where vortices are created and observed 20] . 

In summary, we have examined the vortex core struc- 
ture in the FFLO state in a superconductor under ap- 
plied field. By a simple topological reason, the spin 
polarization of excess up-spin electrons vanishes around 
the place where the vortex line intersects with the nodal 
plane, both of which accommodate excess up-spin elec- 
trons. This intricated spatial magnetization profile may 
be probed by neutron or X-ray diffraction experiments 
as extra-diffraction spots, or NMR and /iSR experiments 



through the resonance line shape. This topological ex- 
planation is confirmed by detailed computations based 
on the BdG equation. We provide full account of these 
features through analyses of the LDOS at various places. 
These spectral features should be directly checked by 
STM-STS experiments. We have proposed these exper- 
iments in connection with the "novel" high field phases 
reported in CeCoIns. The present conclusion is also hold 
for ultra-cold atomic Fermion condensates with unequal 
two species populations under rotation Q. 

We thank H. Adachi, Y. Matsuda, E. M. Forgan and 
H. Yoshizawa for useful discussions. 
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